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Modem, state of the art nanomechanical devices are capable of creating spatial superpositions that are mas- 
sive enough to begin to experimentally access the quantum to classical crossover, and thus force us to consider 
the possible ways in which the usual quantum dynamics may be affected. One recent theoretical proposal de- 
scribes the crossover from unitary quantum mechanics to classical dynamics as a form of spontaneous symmetry 
breaking. Here, we propose a specific experimental setup capable of identifying the source of unitarity break- 
ing in such a mechanism. The experiment is aimed specifically at clarifying the role played by gravity, and 
distinguishes the resulting dynamics from that suggested by alternative scenarios for the quantum to classical 
crossover We give both a theoretical description of the expected dynamics, and a discussion of the involved 
experimental parameter values and the proposed experimental protocol. 



INTRODUCTION 

Experimentalists are pushing to cool cantilevers with low 
intrinsic damping to ever lower temperatures, both in the con- 
text of magnetic resonance force microscopy (MRFM) H], 
and in the quest to cool a mechanical resonator to its zero point 
motion 12^01 and eventually bring it into a quantum mechan- 
ical superposition llj-El. One such experiment has recently 
succeeded in cooling a resonator to its quantum mechanical 
ground state and has demonstrated the controllable creation of 
single quantum excitations [2\. Several other feasible experi- 
ments have been proposed in which a mechanical resonator is 
coupled to a quantum system in such a way that, with some 
experimental progress, a superposition of spatially separated 
states of the resonator may be detectable |6, 7|. With the ad- 
vent of such experiments, involving both quantum mechanics 
and mesoscopic objects, the problem of reconciling quantum 
mechanics with classical physics is no longer a purely theoret- 
ical endeavor, but rather becomes an experimental necessity. 

Many theoretical proposals exist for how to explain the ap- 
parent absence of (spatial) quantum superpositons of macro- 
scopic objects in the classical world. The most well known of 
these include the idea that environment induced decoherence 
hides the superposed states from our view |8|; the idea that 
our personal participation in superpositions implies an inabil- 
ity to observe alternative branches of the superposed state of 
the universe |9 1; and the idea that corrections to Schrodinger's 
equation become important at a scale intermediate between 
that of microscopic particles and macroscopic objects ITOl - 
[131 • Recently, it has been pointed out by several authors that 
gravity may play an important role in scenarios of the latter 
kind llinfT4HT6l . One particular very recent proposal posits 
that the transition between quantum dynamics at the micro- 
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scopic scale and classical physics at the macroscopic scale 
may be described as a form of spontaneous symmetry break- 
ing, in direct analogy to the breaking of translational sym- 
metry in macroscopic crystals, rotational symmetry in macro- 
scopic magnets, and so on 1 17 1. The particular symmetry be- 
ing broken in this case is the unitarity of quantum mechani- 
cal time evolution, and the necessary symmetry breaking field 
may be supplied by the subtle influence of gravity at meso- 
scopic length scales. 

Here, we present a specific experiment which we believe is 
capable of producing a superposition that is massive enough 
to test the theory of spontaneously broken unitarity, while al- 
lowing enough experimental control to also differentiate its 
predictions from those of alternative scenarios. The suggested 
setup is based on realistic estimates of experimental parame- 
ters which may be obtained in state of the art MRFM experi- 
ments. The classical object being forced into a superposition 
is the micromechanical resonator which forms the detector 
arm of a typical MRFM experiment. This setup differs from 
similar proposals in the literature 16] |7| because of the entan- 
glement of the resonator with a nearby microscopic spin state, 
which allows a definitive differentiation between the effects of 
decoherence, spontaneously broken unitarity, and alternative 
scenarios, using the specific experimental protocol described 
here. 

In the following, we first describe the proposed experimen- 
tal setup, and discuss the parameters that will have to be ful- 
filled in order for the experiment to be successful. We then 
present numerical simulations of the expected time evolution 
of this setup in the context of spontaneously broken unitarity, 
which explicitly show how gravity may influence the quantum 
dynamics of the resonator. We also show how a pointer basis 
(which defines the possible outcomes of a measurement on a 
quantum system) and Born's rule (which deals with the prob- 
ability that such a measurement yields a given result), are au- 
tomatically recovered from the time evolution of the resonator 
in this scenario. Finally, we turn to a detailed description of 
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the experimental protocol required to distinguish the specific 
time evolution discussed here from both more common distur- 
bances to quantum dynamics, such as decoherence, and from 
alternative theoretical scenarios. 



THE EXPERIMENTAL SETUP 

Inspired by the experiment by Rugar et al. fTSl, in which 
the force exerted by an electron spin on a small magnet is 
detectable by measuring the deflection of a mechanical res- 
onator, we have proposed that such an experimental configu- 
ration can be used to bring a significant mass in a superpo- 
sition involving a large displacement llT6l . The adaptation of 
such a setup shown in figure[T]consists of a thin wire holding 
a plate of mass m as well as a small spherical magnet with 
magnetization M, in close vicinity to a single, isolated elec- 
tron spin. A natural candidate for the electron spin is the well 
known Nitrogen- Vacancy (NV) color centre in diamond be- 
cause already at room temperature it has shown a coherence 
time as large as 250 microsec |19|. In isotopically engineered 
diamond this increases to 1.8 msec [201, while at cryogenic 
temperatures the T2 time may be expected to become even 
larger Rabl et al. have developed a purely quantum mechan- 
ical description of this mechanical resonator coupled to a sin- 
gle electron spin II2TI . For independent detection of both the 
cantilever and the electron spin we envision that the cantilever 
motion can be observed by coupling it to a SQUID through 
a coil, while the electron spin can be characterized optically. 
When no current is injected into the SQUID and when no light 
is coupled to the NV centre, the purely quantum mechanical 
description of Rabl et al. should normafly be apphcable. 

While the non-magnetic = 0) ground state of the spin 
leaves the resonator untouched, the = — 1) (and + 1)) 
excited state will attract (and repel) the cantilever with a force 
F = jxG where G is the magnetic field gradient originating 
from the magnetic sphere at the position of the electron spin 
with magnetization pL. Inducing transitions between and su- 
perpositions of the different spin states can readily be achieved 
by applying an appropriate combination of static Bq and radio 
frequency B\ magnetic fields. Starting from a superposed spin 
state, the magnetic coupling between the spin state and the 
motion of the resonator then eventually yields a superposition 
of two out of phase oscillation modes of the resonator. 

For the experimental protocol described in detail in section 
[j and employing the realistic parameters discussed presently, 
the amplitude of the two anti-phase oscillations could be 
caused to approach the thickness of the plate. Importantly, 
these amplitudes could in principle be achieved in a time 
shorter than the expected dephasing time of the electron spin 
due to nuclear spins in the diamond lattice |20| and shorter 
than the dephasing time of the cantilever due to coupling to the 
phonon bath ll22l|23l . Thus, a superposition of the plate over 
macroscopically distinct positions may be achieved within the 
limits set by decoherence. 




FIG. 1: The proposed experimental configuration. A thin gold plate 
is attached to a nanowire or multi-walled carbon nanotube to create 
a massive mechanical resonator. A small spherical magnet acts as a 
detector for the spin magnetic moment in a diamond NV centre. The 
graph on the right shows the energy levels of the spin as a function 
of magnetic field. 

Order of Magnitude Estimates for Experimental Parameters 

The experimental parameters required to be able to un- 
ambiguously assign an observed decay time to any process 
other than decoherence lie only just beyond what is already 
being used in present day MRFM measurements. In the 
experiments of Rugar et al. the deflection of a resonator 
due to the interaction with a single electron spin was mea- 
sured [18]. The field gradient in these experiments was 
dB/dx = 2-10^ T/m which, with a magnetic moment of 
the electron spin of jig = 9.3 • 10^^^ J/T, leads to a force 
■P'spin = liBdB/dx = 1 .8 • 10"'^ N. With the stiffness of the res- 
onator A: = 1 . 1 • 10^''^ N/m, this would imply a static deflection 
of only c/q ~ Fspin/k — 1.7 ■ 10^'^ m. 

To increase the deflection of the resonator, the electron spin 
was inverted twice during each resonator period, TJ-es. If the 
spin inversions remain coherent with the resonator motion 
long enough, the amplitude could in principle be increased 
to d = QFf^pin/k:, where Q is the Q-f actor of the resonator 
Since in practice the spin inversions remain coherent with the 
cantilever motion only for a time T„, ^ QTies, the maximum 
amplitude is limited to li = (1 — exp(— T,„/27ies)) ■ QFspm/k, 
where T,„ is the so-called rotating frame relaxation time. For 
the experiment by Rugar et al. T,„ = 760 msec and the maxi- 
mum amplitude would exceed 60 pm 1 1 1 . 

For a spin in a quantum superposition of states, T„, would 
have to be replaced by a dephasing time T2, which has been 
measured to be 1 .2 msec in NV centers in diamond in which 
the '^C isotopic content was reduced 1.20.1. Although this is 
the longest decoherence time measured in a solid state system 
to date, the fact that T2 is so much shorter than T„, still limits 
the distance between the two centers of mass involved in the 
superposition. Fortunately, in the context of nuclear magnetic 
resonance force microscopy, the group of Rugar now employs 
field gradients of 4 • 10^ T/m 1 1 1, which when applied to an 
electron spin experiment, would lead to a 20 fold increase of 
the deflections. An alternative way of enhancing the deflec- 
tion would be to employ softer springs. In our lab L24i . we 
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have been able to fabricate nanowire resonators with spring 
constants down to k = I ■ 10^^ N/m. Such a nanowire, ending 
in a thin gold foil measuring 5 jJ-m x 5 jUm x 20 run, would 
have a resonance frequency of 1 .6 kHz. 

Combining an increase of the NV-centre dephasing time to 
T2 = 50 msec (which might be achievable by further elim- 
inating nuclear spins in the diamond or by working at suf- 
ficiently low temperatures to freeze out the nuclear spins) 
with the highest field gradients in the Rugar lab and a soft 
k — 1 ■ 10^^ N/m resonator would result in a deflection of 
d = (72/7'res) • [pLg/k) ■ dB/dx — 3 nm. As explained in the 
next section, the typical energy scale involved in the pro- 
cess of spontaneous unitarity breaking may be expected to be 
A = Gir^-if- jl} = 4 • 10^^^ J (with L the sample thickness and 
G the gravitational constant), which yields a typical time scale 
h/A = 25 msec. Since this is of the same order of magnitude 
as the dephasing time of the NV-centre, it may be possible for 
the loss of unitarity to become apparent before the NV centre 
decoheres. 

The dephasing of the cantilever itself, on the other hand, 
may be an issue. This time scale can be estimated by describ- 
ing the environment as an infinite bath of harmonic oscilla- 
tors and integrating out the environmental degrees of freedom 
II25] |26l . For the specific case where the dephasing is due 
to the phonons which are excited within the resonator, this 
problem has received some attention in the context of grav- 
itational wave detectors and the timescale resulting after the 
phonon bath has been integrated out was shown to be Tres = 
grres/(27rA^phon), where A^phon = kBTyy^ode/hfres is the number 
of phonons of the resonator mode, which is determined solely 
by the resonator's thermodynamic mode temperature Ty^ode, 
and frequency /res- In order for this dephasing time to exceed 
h/A, the temperature of a resonator with resonance frequency 
/res and Q — 10^ would have to be smaller than 100 jxK, 
which would require the experiment to be done in a nu- 
clear demagnetization refrigerator However, although before 
we considered the damping J ^kj {2nfQ) = 1 • 10^'^ Ns/m 
with g = 10^, it is possible that future resonators will have 
even higher Q since much lower damping values have been 
achieved in the carbon nanotube resonators fabricated by Hut- 
tel et al. 127]|, who estimate 7 ^ 27r/reswVfi = 2 • 10"'^. If 
such a low damping can be achieved for the resonator pro- 
posed here, the mode temperature may be around 1 mK with- 
out dominating the dephasing time of the NV-centre. 

It has been argued before that the thermal contribution to 
the dephasing time considered above should give a correct or- 
der of magnitude estimate even in the presence of a driving 
force |6|. However, if we consider a dephasing mechanism 
whose strength depends explicitly on the amplitude of the os- 
cillation as well as its mode temperature [,28,], we should re- 
place A^phon in the expression for Tres by jkd'^/hfres- With the 
lowest damping considered (y — 2- 10^'^) at a temperature of 
100 IJ.K, this results in a dephasing time of the order of 1 msec. 
In that case, one would need to resort to using an even thinner 
gold plate to see the effect of the spontaneously broken uni- 
tarity. A gold mass of 10 nm thickness (f.e. deposited on a 



single sheet of graphene), in the presence of such low damp- 
ing, gives rise to the timescale h/A = 20 msec, which might 
just be detectable within the time limit set by decoherence. 

Although the phonon contribution considered here is ex- 
pected to dominate the dephasing of the resonator, other 
mechanisms are also present. Defects in the cantilever for ex- 
ample may act as effective two-level systems at low temper- 
atures, and the clamping of the cantilever may lead to addi- 
tional phonon radiation, both contributing to the dephasing of 
the resonator |29|. The former issue can be addressed by the 
application of small additional magnetic fields or by reducing 
the number of defects in the carbon nanotube resonator, while 
a scheme to circumvent the problem of clamping loss using 
optical trapping of a resonator has recently been proposed in 
a different context ||30ll . 

SIMULATING THE TIME EVOLUTION 

In this section, we will focus on spontaneously broken uni- 
tarity as one particular scenario for the transition from quan- 
tum dynamics to classical physics, and show explicitly how 
its predictions may be expected to influence the dynamics of 
the coupled resonator-spin system. How the proposed experi- 
mental setup may be used to also distinguish these particular 
predictions from those of other scenarios will be discussed in 
the next section. A detailed discussion of the theory of spon- 
taneously broken unitarity can be found in ref. ifTTll . and we 
will only repeat its main results here. 

The basic observation underlying the idea of spontaneous 
unitarity breaking is the fact that the dynamics generated by 
the usual time dependent Schrodinger equation is unstable, in 
the sense that even an infinitesimally weak perturbation may 
qualitatively change its behavior in the thermodynamic limit. 
This situation is analogous to that of other spontaneously bro- 
ken symmetries: the Hamiltonian for a crystal is invariant un- 
der translations, but the sensitivity to even infinitesimal pertur- 
bations allows a macroscopic crystal to nonetheless localize 
in only a single position (31]. In the same way the unitar- 
ity of quantum mechanical time evolution is sensitive to even 
infinitesimally small non-unitary perturbations, which allows 
sufficiently large objects to undergo classical dynamics |17|. 
An important caveat in this argument is that there has to ex- 
ist a fundamental non-unitary interaction somewhere in nature 
II32II . It has been pointed by several authors that gravity may 
fill this role, and that in fact the energy scale on which the ef- 
fects of gravity as a non-unitary influence become important 
lies in the regime in which we expect the quantum to classical 
crossover to take place |[14l - (T6l . 

Because of the relatively small mass and low density of the 
resonator in the proposed experiment, its dynamics will be 
very close to purely quantum mechanical, and we expect to be 
able to include any non-unitary effects due to gravity as minor 
perturbations to Schrodinger's equation. Following the pro- 
cedure outlined in the appendix, we thus write the dynamics 
of the resonator in the presence of a non-unitary gravitational 



4 



perturbation as: 

where H is the usual quantum mechanical Hamiltonian, G 
is the gravitational constant, and m and L are the mass and 
width of the massive plate. The operator x is the position op- 
erator for the centre of mass of the plate (measured with re- 
spect to the centre of mass of the initial wave function). The 
time-dependent, randomly fluctuating variable £, has been in- 
troduced as a correction to i, because the theory of gravity 
(i.e. general relativity) insists that any quantity which can be 
used as a measure of distance between locations in different 
components of a spacetime superposition must be ill-defined. 

Notice that both general relativity and quantum mechan- 
ics are in fact fully deterministic theories. There is thus no 
reason to believe that any part of the interplay between quan- 
tum mechanics and gravity should be anything but determin- 
istic. The introduction of a random variable in equation ([T]i 
should be seen only as a poor man's approach towards sim- 
ulating an essentially ill-defined quantity: close to the point 
where unitary quantum mechanics is still a good description 
of nature we may get away with using the concept of super- 
positions -even though they really are ill-defined notions in 
general relativity- if we include also an effectively random 
correction to the differences in distance between superposed 
spacetimes. Although superpositions and random variables 
may not actually feature in the exact reality of quantum grav- 
ity, we assume that if we insist on the possibility of effectively 
describing the state of the system as a superposition in some 
limit, then we also need to take into account an effectively ran- 
dom correction to the notion of position. We are thus led to 
the effective, phenomenological description of the first order 
correction to the unitary Schrodinger equation given by equa- 
tion ([T}. In the appendix we give a more detailed discussion 
of both the steps leading to equation ([T]) and the implications 
of the imaginary term on the conservation of energy and on 
the conservation of the norm of the wavefunction. 

In contrast to the Schrodinger-Newton equation, which has 
been formulated in a similar context by other authors [33l,'34|, 
equation ([T} is a purely linear equation, which arises natu- 
rally from the extension of the well studied equilibrium mech- 
anism of spontaneous symmetry breaking to the dynamical 
realm, and obeys all the requirements of a theory of sponta- 
neous symmetry breaking: for unitarity to be spontaneously 
broken, one necessarily needs to invoke a singular thermody- 
namic limit, a macroscopic order paramter, and a symmetry 
breaking field |17|. It is because of this connection to spon- 
taneous symmetry breaking, that the timescale Gm^(x^) /2L^, 
which has been suggested before to set the appropriate time 
scale for the influence of gravity on quantum mechanical time 
evolution [,14-r 16|, may be embedded into the explicit expres- 
sion of the wavefunction dynamics provided by equation ([TJ. 
This final form of the modified Schrodinger equation can be 
straightforwardly integrated using standard numerical meth- 




FIG. 2: The two reduction processes induced by the modified 
Schrodinger equation. Starting at f = with the superposed wave- 
function (x^^ii/'^ {x + d) +Xl¥'^i^-d)) and ^ > 0, the time 
evolution induced by equation l[TJ with a static ^ initially causes the 
fast reduction of the wavefunction to just the component Xi v" {x — d) 
(solid curves), and then the slow reduction of that component towards 
a wavefunction centered at x = ^ (dashed curves). For clarity, all 
curves have been renormalized before plotting them. The parameters 
used in this simulation are based on the optimistic order of magnitude 
estimates discussed in section| They are m = 1- 10"''*kg, L = 1 nm 
and d = I nm, resulting in an initial reduction time of the order of 
t = I ms. 



ods to yield a phenomenological prediction for the expected 
time evolution of our resonator experiment. 



Two Time Scales 

To be explicit, we first express the state of the system 
depicted in figure [T] by the quantum numbers for the posi- 
tion of the cantilever's centre of mass x and the orientation 
of the spin's magnetic moment a. We are then interested 
in the dynamics of the system starting from the initial state 
(p{x,t = 0) = ax^w'^ix + d) + PXiV'^ix — d), where ^^{x) 
is the groundstate wavefunction of the harmonic oscillator of 
mass m centered at x = 0, and Xa indicates the spin state with 
5~ = C7. The entangled initial state is prepared by performing 
an MRFM measurement on a suitably initialized superposi- 
tion of the spin state, as discussed in the next section. 

For the purpose of clarity we ignore the imposed oscillatory 
motion of the cantilever in this simulation and assume instead 
that it is fixed at its maximum displacement ±d. It is straight- 
forward to also include these oscillations but they shroud the 
role of the unusual dynamics induced by the phenomenologi- 
cal gravitational term and do not change our conclusions. Fur- 
ther ignoring for the moment the phonons and other internal 
degrees of freedom of the oscillator, we are left with only the 
kinetic energy of the oscillator as a whole to set the Hamil- 
tonian in equation ([T]|, so that we have H ~ j(2m). The 
spin in the diamond NV centre is assumed to be free for the 
duration of the experiment. 

The time-dependent wavefunction of the resonator found 
by solving the differential equation Q, now shows two dis- 
tinct processes which happen simultaneously. These can be 
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FIG. 3: Two typical realizations of the dynamics of equation |T]( in 
the presence of a randomly fluctuating ^(t), starting from the initial 
wavefunction X-\Mf^{x + d) + \/3/4 Xiv'^i^-d). The fate of 

the wavefunction is tracked by the expectation value of the position 
operator, (x). The initial wavefunction is reduced to just one of its 
components in a time proportional to the inverse of the strength of 
the gravitational term in equation H i. The parameters used in these 
simulations are the same as in figure 2] 
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FIG. 4: The cumulative proportion of wavefunctions being reduced 
to Xi^\i/^{x + d), starting from a X-^W^H-^ + d) + jS XlV^\-^^d). The 
initial dynamics shown in figure[3]was repeated 20000 times for each 
of three different sets of the parameters m and d. For each set the cu- 
mulative proportion of outcomes resulting in XfV^i^ + d) is seen to 
be lap, in accordance with Bom's rule. Here we used a = ^9/20. 



Bern's Rule 



most easily understood for the special case of constant ^ (t). In 
that case the gravitational term of equation ([T]i exponentially 
suppresses the weights of both components of the wavefunc- 
tion, but one component is suppressed more strongly than the 
other, depending on whether x — ^ is positive or negative. One 
component thus quickly dominates the overall wavefunction, 
and the superposition is reduced to just the single product state 
X^yf^ix + d) or XiY^ix — d), as is shown in figure |2] That 
the reduction of the superposed wavefunction results in just 
one component in the position basis is a direct consequence 
of gravity providing the unitarity breaking field in equation 
Q. Because of the special role played by the position variable 
in the theory of gravity, spatial superpositions become unsta- 
ble and the pointer basis which defines the possible outcomes 
of the reduction process consists of spatially localized states. 
Notice that the spatial reduction of the cantilever wavefunc- 
tion automatically implies that also the entangled spin state 
will end up in only one of its components, even though it is 
not itself subject to any gravitational effects. 

The localization dynamics is taken even one step further 
in the second process induced by the modified Schrodinger 
equation. Because the single component \i/^{x±d) is itself 
a wavefunction spread out in space, it too is subject to the 
uneven suppression by the gravitational term. As long as the 
spread of the wavefunction (i.e. its deBroglie wavelength) is 
short compared to the distance £,±d however, the difference 
in amplification rates of its components is very small, and the 
secondary reduction of the wavefunction into a state centered 
aix — ^ will be very slow (see figure[2|i. 



In the presence of a fluctuating ^{t), the same two processes 
occur This time however x — ^ may alternate in time between 
being positive or negative. It will thus alternately amplify the 
relative weights of the two different components of the ini- 
tial superposition, and the question of which component wins 
out in the end becomes a realization of the "gambler's ruin" 
game: each component will randomly increase and decrease 
in relative weight until one of them completely dominates. Al- 
though the fluctuations of the stochastic variable B,{t) slow 
down the rise to dominance of either of the components, a fi- 
nal resolution can still be seen to be reached within a typical 
timescale proportional to the inverse of the gravitational en- 
ergy A = GnP'd^ /L^ defined in equation ([T]) 1 17 1, as expected 
on general dimensional grounds |14-16|. This energy scale 
sets the minimum time scale required by the resonator dynam- 
ics to unambiguously display non-unitary effects. 

An example of a typical solution of the modified 
Schrodinger equation with a fluctuating stochastic variable is 
shown in figure |3] Because ^(f) is a randomly oscillating 
variable, it is impossible to predict which of the two com- 
ponents will survive the imposed dynamics. If one of the 
components has a larger weight in the initial wavefunction 
however, it is more likely to survive. In fact it is possible 
to mathematically prove that after an infinitely long time the 
only possible average outcome of this process is the emer- 
gence of Born's Rule 1321 l35l . That is, if we repeat the sim- 
ulation many times always starting from the same initial state 
aX|i/*'(x + c/) + j3x|V^^(x — c/), but randomly generating a 
new ^(f) for each run, the average fraction of evolutions re- 
sulting in the final state XA^'yfi{x + d) converges to |a|^, as in 
figure]?] 
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PROTOCOL FOR THE EXPERIMENT 

Returning now to the experimental setup described in sec- 
tion! need to find a way of distinguishing the dynamics of 
spontaneous unitarity breaking as described by equation ([T]i 
from both the predictions of alternative models of the quan- 
tum to classical crossover, and from the ubiquitous decoher- 
ence effects. Especially the latter is not an easy task, because 
in the proposed setup, it is not possible to directly observe 
the time evolution of the state of the cantilever as displayed 
in figure [3] and we have to rely on ensemble measurements 
instead. 



Distinguishing Decoherence from Decay 

The central observation in the theory of decoherence is 
that the many unobservable microscopic degrees of freedom 
which are necessarily present in any experiment, have the ef- 
fect, if averaged over an ensemble of measurements, to re- 
duce any initially pure density matrix to a classical mixture 
of pointer states [8|. Although decoherence and the non- 
unitary dynamics of equation ([T]i have completely different 
origins and differ substantially in their descriptions of a sin- 
gle measurement, the ensemble averages look very similar 
In terms of the reduced density matrices (i.e. density matri- 
ces averaged over all unobservable environmental degrees of 
freedom), both give rise to an effective, ensemble-averaged 
evolution of the type 
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Here T is either the decoherence time Tdecoh or the average 
decay time Tdecay To distinguish the influence of the non- 
unitary term from the more usual effects coming from various 
sources of decoherence, we propose an experimental protocol 
consisting of two different stages. 



Stage One 

In the first instance, the spin and the resonator should be 
entangled in such a way as to give rise to a spatial superpo- 
sition of resonator states. This can be achieved in a standard 
Magnetic Resonant Force Microscopy experiment in which 
the spin is reversed at twice the resonator frequency to drive 
the resonator, but starting from a superposed spin state. One 
component of the spin wavefunction will then enhance the res- 
onator motion more strongly than the other, and as long as 
coherence can be maintained, the result will be an entangled 
state in which the resonator motion has a larger amplitude in 
one component than in the other: 

(p{t = 0) = axt¥^{x) + liXi¥H^) 



^p(f = 0) 



a 



where V^^(x) and V'^(x) indicate the wavefunctions of the os- 
cillator with enhanced or suppressed amplitude respectively, 
and the spin's S' component in Xa is measured in the coro- 
tating frame. This state with spatially separated components 
for the massive resonator will be sensitive to gravity-induced 
non-unitary effects as well as to decoherence, so that the off- 
diagonal element of its reduced density matrix after a given 
time interval ti is given by: 



Pl2(fl) 



(4) 



(3) 



The effective dephasing time (which combines the gravita- 
tional decay time and the effects of decoherence) can be found 
after many repetitions of the experiment in which the final 
state is read out by either coupling the nanomagnet attached 
to the oscillator to a SQUID through a coil, or alternatively by 
measuring the fluorescence coming from the spin at the NV 
centre after the application of a microwave % /2 pulse (which 
rotates the spin orientation into the xy plane). The effective 
dephasing time found in this way, in general depends on both 
the mass and the shape of the resonator The characteristics of 
this dependence can be used as a first means of distinguishing 
gravity-induced non-unitary dynamics from any given, known 
source of decoherence. For any given source of decoherence 
it is straightforward to work out the functional dependence of 
Tdecoh on the mass of the resonator, its linear size, its geomet- 
ric shape, etc. In general, one or more of these dependencies 
will differ from the ones implied by the non-unitary form of 
equation ([TJ. For example, because they depend only on inter- 
nal degrees of freedom of the oscillator, most sources of de- 
coherence are independent of the angle between the direction 
of oscillation and the normal to the surface of the gold plate. 
By measuring the functional dependence of the observed T on 
such parameters, we can thus rule out most known sources of 
decoherence. To make sure that no unknown sources of deco- 
herence are at play either, one can then add a second stage to 
the experiment. 



Stage Two 

Any source of decoherence, no matter what its precise 
physical origin is, will arise from the entanglement of the can- 
tilever motion with some other, unobserved degrees of free- 
dom (called the bath). The gravitational decay on the other 
hand, only requires the presence of a massive superposition, 
and takes place even in complete isolation. To differentiate 
between the two mechanisms, one thus needs to ascertain that 
the cantilever is not entangled with any bath degrees of free- 
dom, known or unknown. This is done in the next state, where 
we again start by creating the entangled superposition of equa- 
tion ([3]l, but rather than reading out the wavefunction of the 
resonator after a given time interval, we imprint its state back 
onto the diamond spin. 

To do this, one could use a spin echo setup, and flip the spin 
at a given instant in time (using an extra microwave n pulse, 
which inverts the orientation of the spin), so that from then 
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on the magnetic force between spin and resonator will tend to 
damp its motion if it enhanced it before the spin flip, and vice 
versa. Thus the effect of the different spin orientations on the 
oscillator position will be effectively undone. Once the res- 
onator has traced its way back to its original wavefunction, the 
position of the resonator and the spin direction are no longer 
entangled and the magnetic driving field can be switched off. 
In the absence of any decoherence and decay processes, the 
evolution of the wavefunction, including the creation of the 
entangled state, would be simply: 

(p(f-0)= [axt + PXil^ix) 
(p{0<t<ti) = axtW^ix,t) + l5XiW^{x,t) 
(pih <t<t2) = aXiW^ix,t) + l5xtW^{x,t) 

(pit = f2) = [axi + Pxt] (5) 

Here the spin echo pulse was applied at f = f i , and the time 
dependence of y/^-^{x,t) signifies a growing (or shrinking) os- 
cillation amplitude as a function of time for t <t\, but the op- 
posite for t > t\, due to the opposing influence of the flipped 
spin. If a non-unitary term is present, the components a and 
j3 become time dependent, but the oscillator state will still be 
brought back to ^^{x) at f = f2. Likewise, decoherence may 
result from the entanglement of these states with bath degrees 
of freedom, but these do not alter the final oscillator state. 

After decoupling the spin state and the oscillator state at f2, 
there will not be any additional gravitational decay because 
there is no longer any spatially superposed mass. Any de- 
coherence caused by the entanglement of the diamond spin 
with unknown internal (bath) degrees of freedom of the res- 
onator however, will continue unabated even if the resonator 
remains at rest. After all, due to the rigidity of the resonator 
(and the associated Goldstone theorem), its internal degrees of 
freedom are decoupled from its collective coordinates. Even 
though the spin flip at t = ty effectively reverses the time 
evolution of the collective, centre-of-mass motion of the res- 
onator, the dynamics of the internal, bath degrees of freedom 
is not reversed, and they remain in an entangled state even af- 
ter t = t2- The full wavefunction evolution in the presence of 
both decay and bath degrees of freedom can thus schemati- 
cally be written as: 

(p(f = 0) = [a(0)zt + j3(0)zj /W^^O) 
9(0 < f < fi) ^a{t)x^^^{x,t)<p^{t)+fi{t)xiwHx,t)<^^{t) 
(pih <t<t2) =a{t)XiW^{x,t)<l)^{t) + I5{t)x^w^ix,t)<l)^{t) 

(Pit > t2) = [a{t2)Xi^'^{t) + l5{t2)Xt^'{t)] /W. (6) 

Here ^^"^(f) represent the evolution of the bath degrees of 
freedom, as influenced by the different amplitude states of 
the oscillator. Their time evolution continues beyond t — t2, 
while the components a and j3 become time-independent after 
the centre of mass state has returned to ij/^ix). After tracing 
out the bath degrees of freedom and simultaneously averaging 
over many realizations of the stochastic dynamics imposed by 



the non-unitary time evolution, the off-diagonal element of the 
reduced density matrix for the final state yields: 

Pnit > t2) = a(0)j3(0)*e"'2(^''''''-°h + l/Tdec.y)g-('-f2)/Tdecoh^ 

(7) 

At the end of the second experimental stage then, the co- 
herence of the diamond spin can again be monitored via op- 
tical readings in many repetitions of the experiment. The ini- 
tial suppression of the off-diagonal matrix elements by the 
factor e"'2/'^dccoh-'2/Tdccay known from the first stage of the 
experiment. The remaining suppression of the off-diagonal 
matrix elements by the factor e^('^'2)/i^dccoh can be attributed 
purely to sources of decoherence. It will most likely be dom- 
inated by the well-characterized decoherence of an isolated 
diamond spin due to sources within the diamond environment 
itself. For the experimental parameters proposed here, this 
timescale should be in the millisecond range, and can be eas- 
ily recognized. Any remaining sources of decoherence which 
had not been identified and eliminated in the first experimen- 
tal stage can now be recognized as an apparent contribution to 
g-('-'2)/Tdccoh^ in excess of that expected from the known, pris- 
tine diamond environment. Once these remaining unidentified 
bath degrees of freedom have been eliminated (f.e. through 
additional cooling), the dephasing measured in the first stage 
of the experiment can be conclusively attributed to a decay 
process instead of decoherence. 

Distinguishing Different Decay Models 

After ruling out decoherence as the source for the suppres- 
sion of off-diagonal matrix elements, we still need to ascertain 
whether the observed effective decay time T is due to the pro- 
cess of spontaneous unitarity breaking, which we focussed on 
in this paper, or to some other non-unitary mechanism. The 
most well-known alternatives of this kind fall in the GRW and 
CSL classes lfT0HT2ll . Other proposed scenarios for the quan- 
tum to classical crossover, such as hidden variable and many 
worlds theories are fully unitary ||9]|36l, and the observation 
of any non-zero decay time in the proposed experiment would 
thus suffice to rule out their relevance in the description of the 
resonator dynamics. 

The main assumption in both the CSL and the GRW theo- 
ries is the existence of a specific process beyond the realm of 
applicability of Schrodinger's equation. In the GRW theory 
this process is assumed to be the spontaneous and instanta- 
neous spatial localization of an elementary particle at partic- 
ular intervals in time |10|. If the frequency with which these 
localization events occur is low enough, it will take an im- 
measurably long time (on average) for any individual particle 
to undergo such an event. Within an extended object consist- 
ing of a macroscopic number of particles on the other hand, 
it will not take long before at least one of the particles is lo- 
calized. Assuming that the extended object possesses some 
degree of rigidity, the localization of a single particle within 
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it will suffice to give the entire wavefunction a definite posi- 
tion in space. Macroscopic objects are thus rapidly reduced 
to position eigenstates, while microscopic particles are free to 
spread out in space. Born's rule can be built into this construc- 
tion by making additional assumption about the localization 
process ifTOil . 

Because the localization process proposed by GRW acts on 
individual microscopic particles, and does not depend on the 
overall state of the system, its predicted decay time scales only 
with the number of particles involved in the massive superpo- 
sition. A dependence of the decay time Tdecay in equation (|7]i 
on the involved mass or shape of the resonator will thus clearly 
indicate a departure from the predictions of the GRW theory, 
and require additional non-unitary effects. 

The CSL (or Continuous Spontaneous Localization) mod- 
els can be seen as extensions of the original GRW model in 
which the addition to quantum mechanics is no longer a set of 
instantaneous localization events, but rather a continuous pro- 
cess which constantly acts to gradually localize the individual 
particles. In most modern CSL models, a smeared mass den- 
sity is taken to be the variable which specifies both the rate and 
the final states of the reduction process ilTI - [l3l . The obtained 
average rate of localization ensures that microscopic particles 
take an immeasurably long time to localize, while macro- 
scopic objects again localize almost instantaneously due to 
their rigidity. Born's rule can emerge spontaneously if the 
localization process is stochastic, which is modeled by the in- 
clusion of white noise in the dynamical equations 1 1 1 1. 

The final form of the reduction process described by the 
CSL model resembles the dynamics of spontaneous unitarity 
breaking in equation ([T]i in various ways. Both predict a re- 
duction of the superposed initial state within a timescale pro- 
portional to the square of the total mass of the involved object; 
both involve a random variable which ensures that Born's rule 
can be obeyed; and both mechanisms select the position ba- 
sis as the pointer basis to which macroscopic objects must be 
reduced. 

However, there are also clear differences between the two 
approaches. For example, as discussed in the review articles 
of Bassi and Ghirardi 1 1 1 1, and of Pearle 1 13 1, the effective re- 
duction rate in the CSL model is set directly by the spread of 
the mass density distributions in different components of the 
initial wavefunction, independent of the relative locations of 
their centres of mass (if the distributions are initially well sep- 
arated). Thus, in contrast to equation ([T]) the reduction rate of 
a macroscopic superposition is independent of the relative dis- 
tance between the massive bodies in its components, as long 
as they are initially non-overlapping. Such a difference would 
be obvious in the experiment, if we could construct a super- 
position of the massive plate over distances greater than its 
thickness. If on the other hand the two components of the 
resonator wavefunction do have a substantial spatial overlap 
throughout the experiment, we need to rely on the dependence 
of the decay time on the precise geometry of the massive plate 
to differentiate between the predictions of the CSL model and 
spontaneous unitarity breaking. In the CSL model, the ef- 



fective decay time can be shown to be proportional to d/L^, 
where d is the thickness of the non-overlapping part of the in- 
volved mass distributions |flT|. From equation ([T]i, it can be 
seen that spontaneous unitarity breaking predicts the rate to be 
proportional to d^ jl} instead 1 16|. 

After ruling out decoherence as the cause for the decay of 
the off diagonal matrix element in the two-stage protocol, we 
can thus use the dependence of the observed decay time on 
the mass and geometry of the proposed setup to conclusively 
distinguish the effects of all currently proposed descriptions 
of the quantum to classical crossover 



CONCLUSIONS 

We have proposed a nanoscale experiment, realizable with 
present-day technology, which can be used to produce an en- 
tangled state of a Nitrogen- Vacancy spin in diamond and a 
mechanical resonator Because the entangled state involves 
the superposition of a massive cantilever over a sizable dis- 
tance, we are forced to consider the effects that a quantum to 
classical crossover might have on the dynamics of this state. 
The proposed experiment is focussed on testing one particular 
scenario for this crossover, which describes the disappearance 
of unitary quantum dynamics in the thermodynamic limit as 
a form of spontaneous symmetry breaking, akin to the spon- 
taneous symmetry breaking observed in macroscopic crystals, 
magnets, etc. By constructing an effective, phenomenological 
description of the proposed experiment, it is shown that the 
predicted effects of spontaneous unitarity breaking on the time 
evolution of the resonator state can be experimentally tested. 
A specific experimental protocol is proposed in order to distin- 
guish the effects of spontaneous unitarity breaking from those 
of both decoherence and other proposed models for the quan- 
tum to classical crossover. 

The phenomenological description of the resonator dynam- 
ics predicts the reduction of a superposed initial state to just 
one of its components to take place on a timescale which lies 
just within experimental reach. The reduction of the cantilever 
state also implies the reduction of the entangled spin state, 
even though the spin itself is not directly under the influence 
of the broken unitarity. The emergence of both a pointer basis 
and Bom's rule for the reduction process can be seen to be 
a direct consequence of the influence of gravity on quantum 
mechanics. 



APPENDIX A: THE NORM OF THE WAVEFUNCTION 

To see how the usual probabilistic predictions of quantum 
mechanics may emerge from the modified time evolution of 
equation ([T]l, even though it does not conserve the normaliza- 
tion of the wavefunction ifTTl . consider a particular run of the 
proposed experiment, with the spin of the diamond NV centre 
initially in a superposition of S~ states, while the resonator is 
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in a single position state: 

(p{t <0) = [axt + liXi\w^{x) 



(8) 



Here, without loss of generality, we assume \a\^ + |j3p = 1. 
This state is a stable state, in the sense that the diamond NV 
centre is not affected by the non-unitary term of the modi- 
fied Schrodinger equation since it does not involve a massive 
superposition, while the oscillator is in a state with a single, 
well-defined position, and is thus subject only to the extremely 
slow dynamics indicated by the dashed curves in figure|2] For 
the estimated experimental parameters discussed in section 11, 
the slow motion of the oscillator happens on an unmeasurably 
long timescale. If we were to consider even heavier objects 
(like for example the pointer on a regular measurement ma- 
chine), this time scale becomes even longer. For all practical 
purposes the product state (p{t < 0) is thus a static configura- 
tion, with a conserved norm. 

In the initial phase of the experiment, the spin state and the 
oscillator wavefunction become entangled by the coupling of 
the spin orientation to the deflection of the cantilever in the 
usual MFRM setup. If this can be done within a time that is 
short compared to both the intrinsic coherence time of the spin 
and the typical decay time imposed by the non-unitary field, 
the resulting state will be given by: 



(p{t^0) = axtW^{x + d) + l5xiW%x-d). 



(9) 



This state is an entangled state, which involves the superposi- 
tion of a massive object over different locations. It is therefore 
subject also to the faster non-unitary decay dynamics (indi- 
cated by solid curves in figure [2|i. For the experimental pa- 
rameters proposed in section II, the fast non-unitary dynamics 
takes milliseconds to complete, and is within the measurable 
regime. For heavier objects, this time scale quickly becomes 
unmeasurably short. If we were to consider a truly macro- 
scopic superposition (involving for example the pointer on a 
regular measurement device), the decay would, for all practi- 
cal purposes, be instantaneous. 

During the decay dynamics, the norm of the wavefunction 
is not conserved. The result, if the decay happens on a time 
scale short compared to the involved intrinsic decoherence 
times, will be either one of two possible outcomes: 



^(f >0)= { or 



NX\'W^{x + d) 
NXiW\x-d), 



(10) 



where the norm depends on the precise dynamics, and is 
certainly not equal to one. Despite the lack of a normalized 
wavefunction, there can be no doubt about the physical in- 
terpretation of the final states in equation ( [TO] i. Each of the 
possible outcomes is a single product state which combines 
an eigenfunction of the z-projection operator of the spin in the 
diamond NV centre with a single, well-defined value for the 
spatial location of the oscillator. It must thus be interpreted 
as representing a spin which is fully oriented along the z-axis. 



and an oscillator in a (classical or generalized coherent) state 
centered at position x±d. The position of the oscillator can 
be registered and recorded, and thus effectively serves as a 
measurement of the spin state. Notice that the final state is 
once again a stable configuration which is insensitive to fur- 
ther non-unitary influences. As long as no sufficiently mas- 
sive superpositions are created, its further time evolution will 
be unitary for all practical purposes, and its norm will be 
conserved. 

We can repeat the above experiment many times, and each 
time record which of the two states in equation ( [TO] l is the 
outcome of that particular realization. The number of experi- 
ments resulting in a state with the oscillator centered atx + d 
divided by the total number of conducted experiments will 
then be seen to converge to the value |ap, as indicated in fig- 
ure |4] Since it is impossible to predict which of the two out- 
comes will be realized in any one particular experiment, this 
is equivalent to having a probabilistic reduction of the initial 
state to just one of its components, with the probability for 
obtaining a given component set by its squared weight in the 
initial wavefunction: 



X^\l/°{x + d), P=\a\ 
Xiw'ix-d), P = \I5\ 



(11) 



Notice that for a general initial norm \No\^ = |ocp + |/3p, we 
would find probabilities |a^|/|A^o|^ and \P\'^/\No\'^ instead. 

As noted before, the reduction process would happen un- 
measurably fast if we were to employ any of the truly massive 
devices that are customarily used as quantum measurements 
machines. In that case, the dynamics described here thus 
reproduces the instantaneous 'collapse of the wavefunction' 
that is traditionally postulated as an addition to Schrodinger's 
equation. For devices much lighter than the proposed mi- 
cromechanical oscillator on the other hand, the non-unitary re- 
duction process takes an unmeasurably long time to complete, 
and the dynamics remains effectively unitary, as prescribed 
by Schrodinger's equation. Both the unitary quantum dynam- 
ics of microscopic particles, and the non-unitary behavior of 
classical objects thus emerge naturally from the dynamics of 
equation ([T]). In particular, the usual probabilistic predictions 
of quantum mechanics are recovered, in spite of the absence 
of a conserved norm of the wavefunction. 



APPENDIX B: THE PERTURBATION TO SCHRODINGER'S 
EQUATION 

The dynamical Schrodinger equation can be seen to be 
subject to spontaneous symmetry breaking, in the sense that 
its result may be qualitatively affected in the thermodynamic 
limit by only an infinitesimally weak unitarity breaking field. 
This description of broken unitarity is a straightforward exten- 
sion of the well known mechanism of spontaneous symmetry 
breaking in crystals, magnets, and so on 1 17 1. For such spon- 
taneous symmetry breaking to be effective in large, but finite. 
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objects, there must exist a small, but finite, non-unitary contri- 
bution to its dynamics. In this appendix we give a brief sketch 
of how a non-unitary such as the one in ([TJ may be seen to 
naturally arise from considerations of the interplay between 
quantum mechanics and general relativity 1 17|. 

As it turns out, the fundamental building blocks of these 
theories (unitarity for quantum mechanics and general covari- 
ance for general relativity) are mutually incompatible con- 
cepts. All theories trying to bridge the gap between general 
relativity and quantum mechanics have to abandon either one 
or both of these principles. For example, the background 
dependence of non-equilibrium string theory breaks general 
covariance, while an equilibrium formulation in Euclidean 
spacetime is necessarily dissipative and thus non-unitary. The 
conflict can be made apparent by comparing simple examples 
of time evolution in the two theories 1 15 1. For the purpose of 
this paper we remain agnostic about the ultimate role of gen- 
eral covariance, but we assume that unitarity is not an exact 
symmetry of nature, and may thus be spontaneously broken 
ilTJ. 

The density of the resonator in our experiment is rather 
modest compared to the densities occurring in black holes or 
in high energy physics. The dynamics of the resonator will 
therefore be very close to purely quantum mechanical, and 
we expect to be able to include the effects of general relativity 
as minor perturbations to Schrodinger's equation. Rather than 
identifying the operator d/dt with just the unitary —{i/h)H, 
we therefore introduce a small perturbative term H': 



in — \i/{r,t,a) ^ [H + H'] \l/{r,t,(y), 



(12) 



where the perturbation H' is assumed to be due to the influ- 
ence of general relativity. 

To find the functional form of the perturbation, we start by 
assuming that the defining property of the conflict between 
general relativity and quantum mechanics is the essential non- 
unitarity of the former. Other ingredients (such as a possible 
non-linearity) will be left to higher order terms. The correc- 
tion term H' then cannot be a Hermitian operator, but on gen- 
eral grounds we would like the time evolution to be invertible 
even in the presence of gravity, so that physics makes sense 
if time flows backwards as well as forwards (although it does 
not necessarily have to look the same). One possible way to 
enforce this constraint is by writing 



dt 



(H-iX) w, 



(13) 



where X is a linear and Hermitian operator The fact that the 
time evolution generated in this way does not conserve energy 
(as measured by H) agrees with the lack of a locally conserved 
energy concept in a non-static configuration of general relativ- 
ity. Of course globally, energy should be a conserved quantity, 
and we will have to choose X such that it takes account of 
that restriction. In fact, the presence of an order parameter in 
rigid macroscopic objects can be shown to automatically re- 
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FIG. 5: The massive superposition of a block of width L over a small 
distance x. 



Store global energy conservation in the thermodynamic limit 

El. 

The second step in our argumentation is to relate the 'en- 
ergy' term X to a measure for the extent to which the pres- 
ence of a given quantum superposition is in conflict with the 
requirements of general relativity, and thus gives rise to a non- 
unitary correction to quantum mechanics. It has been shown 
by Diosi lfT4l . and independently by Penrose ifTSl . that there 
exists a covariant way of constructing such a measure. For 
an equal superposition of two different mass distributions pi 
and P2 in the Newtonian Umit, they found this measure to be 

miiia: 



-47tG 



[piW-p2W][Pi(y)- 



-y\ 



(14) 



In the special case in which the superposed mass m is a block 
which is evenly superposed over a distance x small compared 
to its width L (as in figure |5]l, the integrals of equation ([14]) 
can be evaluated to yield tl6l 



(15) 



„ in 

A = G—.x^ 
2L3 



The form of equation ( fTS) allows a straightforward general- 
ization to the case of a generic superposition which consists of 
any number of components carrying arbitrary weights in the 
wavefunction. If i is the standard quantum mechanical opera- 
tor measuring the position of the block's centre of mass (with 
the zero of position at the overall centre of mass of the ini- 
tial wavefunction), then we can interpret A as the expectation 
value of the quantum operator 



(16) 



Notice that although we started with a semi-classical defini- 
tion of the incompatibility measure in equation ([T4]i, the final 
expression of equation ([T6| is a fully quantum mechanical, 
Hermitian operator. This quantized form is only applicable 
to the specific case of a massive block superposed over small 
distances along one spatial direction, and in that sense is not 
as general as the original semiclassical expression. However, 
because it is an operator form it can be applied to any such 
wavefunction of the massive block, and it is not restricted to 
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an even distribution over just two states. For such a more gen- 
eral wavefunction, the expectation value of the operator is still 
a good measure of the uncertainty introduced into the concept 
of a locally conserved energy by the application of general 
covariance. These properties of the operator form of A sug- 
gest that it may also be an appropriate form for the first order 
correction to Schrodinger's time evolution in equation ( [T3| ): 

2 

Xc<G^x\ (17) 

Finally, we need to address the fact that the spatial distance 
between points in separate space-times is an ill-defined no- 
tion in general relativity ifTSl . In this regard, the expression 
of equation ([TTJi for the first order correction is not fully satis- 
factory. The operator x measures the position of the centre of 
mass in a component of the massive superposition using a uni- 
versal coordinate system which is applied equally to all other 
components. We needed to define such a 'best possible' match 
between the coordinate systems of different components to 
first calculate the semi-classical measure for its inappropriate- 
ness (equation ([T4]l), but the final form of the correction to 
Schrodinger's equation should not depend on any such (ar- 
bitrary) specific choice of coordinate system. To model the 
ill-definedness in the definition of x we introduce a random 
variable £, , and replace the measure of distance by (i — ^ ) . 

X = G^[x-^f (18) 

As emphasized also in the main text, the introduction of a 
random variable here should be seen only as a poor man's 
approach towards modeling an essentially ill-defined quantity. 
We do not expect the full theory of quantum gravity to be 
non-deterministic, but we are forced to take into account an 
effectively random correction to the notion of position if we 
are to consider gravity's first order perturbation to quantum 
dynamics. The equality of the numerical pre-factors of i and 
£, reflects the fact that the correction due to gravity should 
depend on the spatial spread of the wavefunction. 

We thus finally reproduce equation ([T]i as the effective, phe- 
nomenological description of the minimal correction to the 
unitary Schrodinger equation due to the influence of general 
covariance: 

1-^-4 '''' 

This final form of the modified Schrodinger equation can be 
straightforwardly integrated using standard numerical meth- 
ods to yield a phenomenological prediction for the expected 
time evolution of the proposed resonator experiment. 
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